home *** CD-ROM | disk | FTP | other *** search
/ Ham Radio 2000 #2 / Ham Radio 2000 - Volume 2.iso / HAMV2 / MISC / HCAL-27 / TRAPDSGN.BAS (.txt) < prev    next >
Encoding:
GW-BASIC  |  1997-01-29  |  14.6 KB  |  456 lines
  1. 10  'TRAPDSGN - 07 SEP 92 - adapted from HARTRPDI, by  A.C.(AL) Buxton W8NX
  2. 20  '                                                  2225 Woodpark Rd.
  3. 30  '                                                  Akron, Ohio  44333
  4. 40  '                                                  Tel. 1-216-836 3854
  5. 50  '............rev. 28 SEP 96
  6. 60  IF EX$=""THEN EX$="EXIT"
  7. 70  CLS
  8. 80  COLOR 7,0,1
  9. 90  PI=3.14159
  10. 100  U$="##,###.###"
  11. 110  UL$=STRING$(80,205)
  12. 120  '
  13. 130  '.....start
  14. 140  CLS
  15. 150  COLOR 15,2,1
  16. 160  PRINT " TRAP DIPOLE - 3 BAND SINGLE TRAP";
  17. 170  PRINT TAB(57);"by A.C.(AL) Buxton W8NX ";
  18. 180  COLOR 1,0:PRINT STRING$(80,223);
  19. 190  COLOR 7,0
  20. 200  PRINT " Press number in < > to:
  21. 210  PRINT UL$;
  22. 220  PRINT "  < 1 >  RUN program"
  23. 230  PRINT "  < 0 >  EXIT"
  24. 240  PRINT UL$;
  25. 250  Z$=INKEY$
  26. 260  IF Z$="0"THEN CLS:RUN EX$
  27. 270  IF Z$="1"THEN 410
  28. 280  GOTO 250
  29. 290  '
  30. 300  '.....L/C trap
  31. 310  CLS
  32. 320  PRINT "In the program to follow, start your design procedure by selecting"
  33. 330  PRINT "the L/C Tuned Circuit option from the menu, then enter only your"
  34. 340  PRINT "lowest operating frequency and a reactance of 180 ohms. Then use"
  35. 350  PRINT "the coil design options to select a coil."
  36. 360  PRINT
  37. 370  PRINT "Press any key to continue.............."
  38. 380  IF INKEY$=""THEN 380
  39. 390  CLS:RUN"\bas\prog\tunecct"
  40. 400  '
  41. 410  '.....start
  42. 420  PRINT " Do you want the program preamble?  (y/n)"
  43. 430  P$=INKEY$
  44. 440  IF P$="Y"OR P$="y"THEN COLOR 7,0,0:GOTO 480
  45. 450  IF P$="N"OR P$="n"THEN COLOR 7,0,0:GOTO 620
  46. 460  GOTO 430
  47. 470  '
  48. 480  '.....preamble
  49. 490  CLS
  50. 500  OPEN"I",1,"\data\docfiles\trapdsgn.txt"
  51. 510  LN=0
  52. 520  IF EOF(1)THEN 590
  53. 530  LN=LN+1
  54. 540  INPUT#1,A$
  55. 550  PRINT "  ";A$;
  56. 560  IF LN<24 THEN PRINT "":GOTO 580
  57. 570  LN=0:GOSUB 4430:CLS   'screen dump
  58. 580  GOTO 520
  59. 590  CLOSE
  60. 600  GOTO 4380    'end
  61. 610  '
  62. 620  '.....input data
  63. 630  CLS
  64. 640  '.....graphics subroutine
  65. 650  LINE INPUT "ENTER: Today's date.......................................? ";D$
  66. 660  IF D$=""THEN LOCATE CSRLIN-1:GOTO 650
  67. 670  INPUT "ENTER: Lowest operating frequency in MHz..................";FL
  68. 680  IF FL=0 THEN LOCATE CSRLIN-1:GOTO 670
  69. 690  INPUT "ENTER: Medium operating frequency in MHz..................";FM
  70. 700  IF FM=0 THEN LOCATE CSRLIN-1:GOTO 690
  71. 710  INPUT "ENTER: Harmonic operating frequency in MHz................";FHAR
  72. 720  IF FHAR=0 THEN LOCATE CSRLIN-1:GOTO 710
  73. 730  INPUT "ENTER: The order of the odd harmonic frequency............";ORDER
  74. 740  IF ORDER=0 THEN LOCATE CSRLIN-1:GOTO 730
  75. 750  INPUT "ENTER: Diameter of wire radiating elements in inches......";DW
  76. 760  IF DW=0 THEN LOCATE CSRLIN-1:GOTO 750
  77. 770  PRINT "PRESS  number in < > to select:"
  78. 780  PRINT "                 <1>  Coaxial traps
  79. 790  PRINT "                 <2>  Open Inductor-Capacitor traps
  80. 800  Z$=INKEY$
  81. 810  IF Z$="1"THEN TT$="C":GOTO 840
  82. 820  IF Z$="2"THEN TT$="O":GOTO 840
  83. 830  GOTO 800
  84. 840  INPUT "ENTER: Estimated outside diameter of trap (inches)........";DT
  85. 850  IF DT=0 THEN LOCATE CSRLIN-1:GOTO 840
  86. 860  INPUT "ENTER: Estimated length of trap (inches)..................";LH
  87. 870  IF LH=0 THEN LOCATE CSRLIN-1:GOTO 860
  88. 880  INPUT "ENTER: Estimated Q of trap................................";QTRAP
  89. 890  IF QTRAP=0 THEN LOCATE CSRLIN-1:GOTO 880
  90. 900  FPSEUDO=FHAR      'harmonic freq. antenna would have if trap were removed
  91. 910  '
  92. 920  '****************************************
  93. 930  '.....start of iterative loop
  94. 940  '****************************************
  95. 950  CONVCHK=0         'restet convergence check flag
  96. 960  FOR AB=1 TO 10:   'design goes through 10 iterations
  97. 970  KA=1              'initial assumed value for dipole length factor is 1
  98. 980  FMHZ=FHAR         'the subroutine iterates to the correct value
  99. 990  '
  100. 1000  GOSUB 4180        'get KA,Z),XCE,CE
  101. 1010                    'KA=  length factor
  102. 1020                    'CE=  equivalent end capacitance
  103. 1030                    'XCE= equivalent end capacitve reactance
  104. 1040                    'ZO=  characteristic impedance of lambda/2 wire section
  105. 1050  ZT=983.57*(ORDER-1+KA)/4/FPSEUDO
  106. 1060                    'ZT=  total length of monopole antenna
  107. 1070  FMHZ=FM           'medium frequency
  108. 1080  '
  109. 1090  GOSUB 4180        'get KA,Z0,XCE,CE
  110. 1100  CS=0.8*SQR(LH*DT)          'stray shunt capacity of trap to free space
  111. 1110  IF TT$="O" THEN CS=1.2*SQR(LH*DT):Z0=1000*SQR(L/C)
  112. 1120                            'LH = trap length
  113. 1130                            'DT = trap diameter in inches
  114. 1140  XCE=-Z0*TAN(KA*PI/2)      'reactance of the end capacitance
  115. 1150  CE=-10^6/2/PI/FM/XCE
  116. 1160  CSP=CS+CE-0.18*SQR(LH*DT)
  117. 1170    'CSP is the combined effect of the trap stray capacity in free space and 1215   'the end capacitance due to fringing of the electrical field. It assumes 1216   'a 9% negative mutual coupling coefficient
  118. 1180  XCSP=-10^6/2/PI/FM/CSP
  119. 1190  KAP=ATN(-XCSP/Z0)*2/PI
  120. 1200    'KAP = total length factor of dipole at medium frequency FM
  121. 1210  Z1=983.57*KAP/4/FM  'length of inner segment
  122. 1220  Z3=ZT-Z1            'length of outer segment
  123. 1230  FMHZ=FL
  124. 1240  '
  125. 1250  '.....convergence check
  126. 1260  IF CONVCHK THEN 1330
  127. 1270  CLS
  128. 1280  COLOR 0,7
  129. 1290  PRINT " CHECKING THE DESIGN CALCULATIONS FOR CONVERGENCE......."
  130. 1300  COLOR 7,0
  131. 1310  PRINT "--------------------------------------------------------"
  132. 1320  CONVCHK=1     'convergence check flag
  133. 1330  '
  134. 1340  GOSUB 4180    'get KA,Z0,XCE,CE
  135. 1350  XCEL=XCE                         'end reactance at frequency FL
  136. 1360  Z0L=Z0                           'chacteristic impedance at frequency FL
  137. 1370  FFUND=983.57/4*0.95/(Z1+Z3)       'fundamental frequency the antenna would
  138. 1380                                     'have if no traps were present
  139. 1390  NORLGTH=2*(Z1+Z3)*FL/983.57      'normalized physical length of the antenna
  140. 1400  RADRES=(31.04*LOG(NORLGTH)+89.72)'radiation resistance of the antenna
  141. 1410                                     'referred to the current maximums
  142. 1420  L=Z3:KV=1:Z0=Z0L:FMHZ=FL:RL=0.000999999:XL=XCEL:A=PI*RADRES/8/Z0L/(Z1+Z3)
  143. 1430  LAMBDA=983.57*KV/FMHZ            'wavelength in feet
  144. 1440  B=2*PI/LAMBDA
  145. 1450  '
  146. 1460  GOSUB 4010   'transmission line calculation
  147. 1470  RINOUTL=RIN:XINOUTL=XIN:XCS=-1E+06/2/PI/FMHZ/CS
  148. 1480  R7=9.9999E-05:X7=XCS*2:R8=RINOUTL:X8=XINOUTL
  149. 1490  '
  150. 1500  GOSUB 4280   'parallel impedance calculation
  151. 1510  RINOUTLP=R9:XINOUTLP=X9
  152. 1520  L=Z1:Z0=Z0L:FMHZ=FL:RL=RADRES/2:XL=9.9999E-05:A=-PI*RADRES/8/Z0L/(Z1+Z3)
  153. 1530  LAMBDA=983.57*KV/FMHZ:B=2*PI/LAMBDA
  154. 1540  '
  155. 1550  GOSUB 4010   'transmission line calculation
  156. 1560  RININL=RIN:XININL=XIN
  157. 1570  R7=9.9999E-05:X7=2*XCS:R8=RININL:X8=XININL
  158. 1580  '
  159. 1590  GOSUB 4280   'parallel impedance calculation
  160. 1600  RININLP=R9:XININLP=X9
  161. 1610  XTRAPLP=-XININLP-XINOUTLP
  162. 1620  LT=XTRAPLP*(1-(FL/FM)^2)/2/PI/FL 'trap inductance UH & CT=trap capacity PF
  163. 1630  CT=1E+06/4/PI^2/FM^2/LT
  164. 1640  FGEN=FHAR:FMHZ=FGEN
  165. 1650  '
  166. 1660  GOSUB 4180     'get KA,Z0,XCE,CE
  167. 1670  NORLGTH=2*(Z1+Z3)*FMHZ/983.57:RADRES=(31.04*LOG(NORLGTH)+89.72)
  168. 1680  RL=9.9999E-05:XL=XCE:L=Z3:A=PI*RADRES/8/Z0L/(Z1+Z3)*1.34
  169. 1690  LAMBDA=983.57/FMHZ*KV:B=2*PI/LAMBDA
  170. 1700  '
  171. 1710  GOSUB 4010     'transmission line calculation
  172. 1720  R7=RIN:X7=XIN:R8=9.9999E-05:X8=10^6/2/PI/FMHZ/(-CS/2)
  173. 1730  '
  174. 1740  GOSUB 4280     'parallel impedance calculation
  175. 1750  XLT=2*PI*FGEN*LT:XCT=-10^6/2/PI/FGEN/CT:RTRAP=2*PI*FGEN*LT/QTRAP
  176. 1760  MTN=SQR((XLT*XCT)^2+(RTRAP*XCT)^2):OTN=ATN(-RTRAP/XLT)
  177. 1770  MTD=SQR(RTRAP^2+(XLT+XCT)^2):OTD=ATN((XLT+XCT)/RTRAP)
  178. 1780  MT=MTN/MTD:OT=OTN-OTD
  179. 1790  RP=MT*COS(OT):XP=MT*SIN(OT)
  180. 1800  RLP=RP+R9:XLP=XP+X9
  181. 1810  R7=RLP:X7=XLP:R8=9.9999E-05:X8=10^6/2/PI/FMHZ/(-CS/2)
  182. 1820  '
  183. 1830  GOSUB 4280    'parallel impedance calculation
  184. 1840  RL=R9:XL=X9:L=Z1
  185. 1850  '
  186. 1860  GOSUB 4010    'transmissopn line calculation
  187. 1870  FPSEUDO=FPSEUDO+0.009*XIN/RIN*FHAR
  188. 1880  PRINT TAB(1);"LT=";USING U$;LT;:PRINT " >H";
  189. 1890  PRINT TAB(21);"CT=";USING U$;CT;:PRINT " pF";
  190. 1900  PRINT TAB(41);"Z1=";USING U$;Z1;:PRINT " ft";
  191. 1910  PRINT TAB(61);"Z3=";USING U$;Z3;:PRINT " ft"
  192. 1920  PRINT
  193. 1930  NEXT AB
  194. 1940  '*********************************************
  195. 1950  'END OF ITERATIVE LOOP
  196. 1960  '*********************************************
  197. 1970  PRINT
  198. 1980  GOSUB 4430     'screen dump option
  199. 1990  '
  200. 2000  CLS
  201. 2010  PRINT "Do you wish impedance calculations at any operating bands of ";
  202. 2020  PRINT "frequencies? (y/n)"
  203. 2030  Z$=INKEY$
  204. 2040  IF Z$="n"OR Z$="N"THEN 3040
  205. 2050  IF Z$="y"OR Z$="Y"THEN 2080
  206. 2060  GOTO 2030
  207. 2070  '
  208. 2080  '.....impedance calculations
  209. 2090  PRINT
  210. 2100  INPUT "ENTER: Center freq.(MHz) of the desired band ";FCEN
  211. 2110  CLS
  212. 2120  PRINT "               TRAP DIPOLE INPUT IMPEDANCE    " :PRINT
  213. 2130  PRINT TAB(3);"FREQ.(MHz)";TAB(20);"INPUT RES.(ohms)";
  214. 2140  PRINT TAB(40);"INPUT REACT.(ohms)"
  215. 2150  PRINT
  216. 2160  FGEN=FCEN:FMHZ=FGEN:NORLGTH=2*(Z1+Z3)*FGEN/983.57
  217. 2170  RADRES=(31.04*LOG(NORLGTH)+89.72):A=PI*RADRES/4/Z0L/(983.57/FGEN/2)
  218. 2180  IF FGEN>1.2*FM THEN 2240
  219. 2190  IF FGEN<0.8*FM THEN 2240
  220. 2200  NORLGTH=2*Z1*FGEN/983.57:RADRES=(31.04*LOG(NORLGTH)+89.72)
  221. 2210  A=PI*RADRES/4/Z0L/(983.57/FGEN/2)
  222. 2220  '
  223. 2230  '.....print loop
  224. 2240  FOR I=1 TO 10
  225. 2250  RADRES=31.04*LOG(NORLGTH)+89.72
  226. 2260  '
  227. 2270  GOSUB 4180     'get KA,Z0,XCE,CE
  228. 2280  RL=9.9999E-05:XL=XCE:L=Z3:FMHZ=FGEN:LAMBDA=983.57/FMHZ*KV:B=2*PI/LAMBDA
  229. 2290  '
  230. 2300  GOSUB 4010     'transmission line calculation
  231. 2310  CS=0.8*SQR(LH*DT):IF TT$="O" THEN CS=1.2*SQR(LH*DT)
  232. 2320  R7=RIN:X7=XIN:R8=9.9999E-05:X8=10^6/2/PI/FMHZ/(-CS/2)
  233. 2330  '
  234. 2340  GOSUB 4280     'parallel impedance calculation
  235. 2350  CS=0.8*SQR(LH*DT):IF TT$="O" THEN CS=1.2*(LH*DT)^0.5
  236. 2360  XLT=2*PI*FGEN*LT:XCT=-10^6/2/PI/FGEN/CT
  237. 2370  RTRAP=0.000999999
  238. 2380  MTN=SQR((XLT*XCT)^2+(RTRAP*XCT)^2):OTN=ATN(-RTRAP/XLT)
  239. 2390  MTD=SQR(RTRAP^2+(XLT+XCT)^2):OTD=ATN((XLT+XCT)/RTRAP)
  240. 2400  MT=MTN/MTD:OT=OTN-OTD
  241. 2410  RP=MT*COS(OT):XP=MT*SIN(OT)
  242. 2420  RLP=RP+R9:XLP=XP+X9
  243. 2430  IF FGEN<0.8*FM THEN 2460
  244. 2440  IF FGEN>1.2*FM THEN 2460
  245. 2450  CS=2*CSP
  246. 2460  R7=RLP:X7=XLP:R8=9.9999E-05:X8=10^6/2/PI/FMHZ/(-CS/2)
  247. 2470  '
  248. 2480  '.....parallel impedance subroutine
  249. 2490  GOSUB 4280
  250. 2500  RL=R9:XL=X9:L=Z1
  251. 2510  '
  252. 2520  GOSUB 4010    'transmission line calculation
  253. 2530  ERRR=(RADRES-2*RIN)/2/RIN
  254. 2540  A=A*(1+0.8*ERRR)
  255. 2550  NEXT I
  256. 2560  '
  257. 2570  FOR J=1 TO 15
  258. 2580  FGEN=FCEN+(J-8)*0.01*FCEN
  259. 2590  FMHZ=FGEN
  260. 2600  '
  261. 2610  GOSUB 4180    'get KA,Z0,XCE,CE
  262. 2620  RL=9.9999E-05:XL=XCE:L=Z3:FMHZ=FGEN:LAMBDA=983.57/FMHZ*KV:B=2*PI/LAMBDA
  263. 2630  '
  264. 2640  GOSUB 4010    'transmission line calculation
  265. 2650  CS=0.8*SQR(LH*DT)
  266. 2660  R7=RIN:X7=XIN:R8=9.9999E-05:X8=10^6/2/PI/FMHZ/(-CS/2)
  267. 2670  '
  268. 2680  GOSUB 4280    'parallel impedance calculation
  269. 2690  CS=0.8*SQR(LH*DT):IF TT$="O" THEN CS=1.2*SQR(LH*DT)
  270. 2700  IF FGEN<0.8*FM THEN 2730
  271. 2710  IF FGEN>1.2*FM THEN 2730
  272. 2720  CS=2*CSP
  273. 2730  XLT=2*PI*FGEN*LT:XCT=-10^6/2/PI/FGEN/CT
  274. 2740  RTRAP=2*PI*FGEN*LT/QTRAP
  275. 2750  MTN=SQR((XLT*XCT)^2+(RTRAP*XCT)^2):OTN=ATN(-RTRAP/XLT)
  276. 2760  MTD=SQR(RTRAP^2+(XLT+XCT)^2):OTD=ATN((XLT+XCT)/RTRAP)
  277. 2770  MT=MTN/MTD:OT=OTN-OTD
  278. 2780  RP=MT*COS(OT):XP=MT*SIN(OT)
  279. 2790  RLP=RP+R9:XLP=XP+X9
  280. 2800  R7=RLP:X7=XLP:R8=9.999E-06:X8=10^6/2/PI/FMHZ/(-CS/2)
  281. 2810  '
  282. 2820  GOSUB 4280   'parallel impedance calculation
  283. 2830  RL=R9:XL=X9:L=Z1
  284. 2840  '
  285. 2850  GOSUB 4010   'transmission line calculation
  286. 2860  '
  287. 2870  '.....screen display
  288. 2880  PRINT TAB(3);USING U$;FGEN;
  289. 2890  PRINT TAB(22);USING U$;2*RIN;
  290. 2900  PRINT TAB(42);USING U$;2*XIN
  291. 2910  NEXT J
  292. 2920  '
  293. 2930  '.....end of impedance calculation loop
  294. 2940  '
  295. 2950  GOSUB 4430     'screen dump option subroutine
  296. 2960  CLS
  297. 2970  PRINT "Do you wish impedance calculations over another band of freqs.?";
  298. 2980  PRINT "  (Y/N)"
  299. 2990  Z$=INKEY$
  300. 3000  IF Z$="Y"OR Z$="y"THEN 2100
  301. 3010  IF Z$="N"OR Z$="n"THEN 3040
  302. 3020  GOTO 2990
  303. 3030  '
  304. 3040  CLS:PRINT TAB(28);"SUMMARY OF SPECIFICATIONS"
  305. 3050  PRINT TAB(28);"  HARMONIC TRAP DIPOLE "
  306. 3060  PRINT TAB((80-LEN(D$))/2);D$:PRINT
  307. 3070  TB=13            'left margin
  308. 3080  PRINT TAB(TB);
  309. 3090  PRINT "Low operating frequency.................. ";USING U$;FL;
  310. 3100  PRINT " Mhz"
  311. 3110  PRINT TAB(TB);
  312. 3120  PRINT "Medium operating frequency................";USING U$;FM;
  313. 3130  PRINT " Mhz"
  314. 3140  PRINT TAB(TB);
  315. 3150  PRINT "Harmonic operating frequency..............";USING U$;FHAR;
  316. 3160  PRINT " Mhz"
  317. 3170  PRINT TAB(TB);
  318. 3180  PRINT "Order of harmonic.........................";USING "######";ORDER
  319. 3190  PRINT TAB(TB);
  320. 3200  PRINT "Length of inner monopole element..........";USING U$;Z1;
  321. 3210  PRINT " ft."
  322. 3220  PRINT TAB(TB);
  323. 3230  PRINT "Length of outer extension beyond trap.....";USING U$;Z3;
  324. 3240  PRINT " ft."
  325. 3250  PRINT TAB(TB);
  326. 3260  PRINT "Overall length of dipole anenna...........";USING U$;2*(Z1+Z3);
  327. 3270  PRINT " ft."
  328. 3280  PRINT TAB(TB);
  329. 3290  PRINT "Diameter of radiating elements............";USING U$;DW;
  330. 3300  PRINT " in."
  331. 3310  PRINT TAB(TB);
  332. 3320  PRINT "Trap resonant frequency...................";USING U$;FM;
  333. 3330  PRINT " MHz."
  334. 3340  PRINT TAB(TB);
  335. 3350  PRINT "Trap inductance...........................";USING U$;LT;
  336. 3360  PRINT " >H"
  337. 3370  PRINT TAB(TB);
  338. 3380  PRINT "Trap capacitance..........................";USING U$;CT;
  339. 3390  PRINT " pF"
  340. 3400  CXR=1000*(LT/CT)^0.5
  341. 3410  PRINT TAB(TB);
  342. 3420  PRINT "Trap characteristic reactance.............";USING U$;CXR;
  343. 3430  PRINT " ohms"
  344. 3440  PRINT TAB(TB);
  345. 3450  PRINT "Estimated outside diameter of trap........";USING U$;DT;
  346. 3460  PRINT " in."
  347. 3470  PRINT TAB(TB);
  348. 3480  PRINT "Estimated length of trap windings.........";USING U$;LH;
  349. 3490  PRINT " in."
  350. 3500  PRINT TAB(TB);
  351. 3510  PRINT "Estimated Q of trap.......................";USING "######";QTRAP
  352. 3520  SCI=0.6*SQR(LH*DT)
  353. 3530  PRINT TAB(TB);
  354. 3540  PRINT "Estimated trap stray capacity to infinite.";USING U$;SCI;
  355. 3550  PRINT " pF"
  356. 3560  L=0.00508*(Z1+Z3)*12*(LOG(4*12*(Z1+Z3)/DW)-1)
  357. 3570  C=7.36*(Z1+Z3)/(LOG(2*12*(Z1+Z3)/DW)/LOG(10)-0.5)
  358. 3580  GOSUB 4430     'screen dump option subroutine
  359. 3590  '
  360. 3600  '.....diagram
  361. 3610  CLS
  362. 3620  PRINT STRING$(79,205)
  363. 3630  LOCATE CSRLIN-1,36:PRINT " W 8 N X "
  364. 3640  PRINT TAB(26);"HARMONIC TRAP DIPOLE ANTENNA"
  365. 3650  PRINT TAB(26);"THENTHENTHENTHENTHENTHENTHENTHENTHENTHENTHENTHENTHENTHENTHENTHENTHENTHENTHENTHENTHENTHENTHENTHENTHENTHENTHENTHEN"
  366. 3660  F$="###.###"
  367. 3670  PRINT TAB(26);"Low operating freq.(MHz)....";USING F$;FL
  368. 3680  PRINT TAB(26);"Mid operating freq.(MHz)....";USING F$;FM
  369. 3690  PRINT TAB(26);"Hi  operating freq.(MHz)....";USING F$;FHAR
  370. 3700  LOCATE  8:PRINT TAB(6);STRING$(70,45)
  371. 3710  LOCATE 8,6:PRINT "<":LOCATE 8,75:PRINT ">"
  372. 3720  LOCATE 8,36:PRINT USING "###.#";2*(Z1 + Z3);:PRINT " ft."
  373. 3730  LOCATE 12:PRINT TAB(6);STRING$(70,196)
  374. 3740  LOCATE 12,39:PRINT "COLOR  VARPTR"
  375. 3750  LOCATE 11,19:PRINT "VARPTRSOUNDORORORSOUNDCOLOR":LOCATE 11,56:PRINT "VARPTRSOUNDORORORSOUNDCOLOR"
  376. 3760  LOCATE 12,19:PRINT "<0xB4!>     BLOAD":LOCATE 12,56:PRINT "<0xB4!>     BLOAD"
  377. 3770  LOCATE 13,19:PRINT "CLSSOUNDSOUNDUSINGSOUNDSOUND'":LOCATE 13,56:PRINT "CLSSOUNDSOUNDUSINGSOUNDSOUND'"
  378. 3780  FOR J=13 TO 17:LOCATE J,39:PRINT "CALL  CALL":NEXT J
  379. 3790  LOCATE 19,33:PRINT "50-72 - Feedline"
  380. 3800  F$="###.#"
  381. 3810  LOCATE 10,19:PRINT USING F$;LT;:PRINT ">H"
  382. 3820  LOCATE 10,56:PRINT USING F$;LT;:PRINT ">H"
  383. 3830  LOCATE 14,19:PRINT USING F$;CT;:PRINT "pF"
  384. 3840  LOCATE 14,56:PRINT USING F$;CT;:PRINT "pF"
  385. 3850  F$="###.##"
  386. 3860  LOCATE 12,11:PRINT USING F$;Z3;:PRINT "'"
  387. 3870  LOCATE 12,27:PRINT USING F$;Z1;:PRINT "'"
  388. 3880  LOCATE 12,48:PRINT USING F$;Z1;:PRINT "'"
  389. 3890  LOCATE 12,64:PRINT USING F$;Z3;:PRINT "'"
  390. 3900  Y=(Z1-INT(Z1))*12
  391. 3910  LOCATE 21,23:PRINT USING F$;Z1;:PRINT "' Inner elements =";
  392. 3920  PRINT USING "###";INT(Z1);:PRINT "'-";USING "##.##";Y;:PRINT CHR$(34)
  393. 3930  Y=(Z3-INT(Z3))*12
  394. 3940  LOCATE 22,23:PRINT USING F$;Z3;:PRINT "' Outer elements =";
  395. 3950  PRINT USING "###";INT(Z3);:PRINT "'-";USING "##.##";Y;:PRINT CHR$(34)
  396. 3960  PRINT STRING$(79,205)
  397. 3970  'GOSUB 4390     'screen dump option subroutine
  398. 3980  GOTO 4380   'end
  399. 3990  END
  400. 4000  '
  401. 4010  '.....transmission line subroutine
  402. 4020  US=(EXP(A*L)-EXP(-A*L))/2*COS(B*L)
  403. 4030  VS=(EXP(A*L)+EXP(-A*L))/2*SIN(B*L)
  404. 4040  UC=(EXP(A*L)+EXP(-A*L))/2*COS(B*L)
  405. 4050  VC=(EXP(A*L)-EXP(-A*L))/2*SIN(B*L)
  406. 4060  C=US-XL*VC/Z0+RL*UC/Z0
  407. 4070  D=VS+RL*VC/Z0+XL*UC/Z0
  408. 4080  E=UC+RL*US/Z0-XL*VS/Z0
  409. 4090  F=VC+RL*VS/Z0+XL*US/Z0
  410. 4100  MNN=SQR(C^2+D^2):THETAN=ATN(D/C)
  411. 4110  IF C<0 THEN THETAN=THETAN+PI
  412. 4120  MND=SQR(E^2+F^2):THETAD=ATN(F/E)
  413. 4130  IF E<0 THEN THETAD=THETAD+PI
  414. 4140  M=MNN/MND:THETA=THETAN-THETAD
  415. 4150  RIN=Z0*COS(THETA)*M:XIN=Z0*SIN(THETA)*M
  416. 4160  RETURN
  417. 4170  '
  418. 4180  '....KA,Z0,XCE,CE subroutine
  419. 4190  FOR W=1 TO 10
  420. 4200  ZX=983.57/FMHZ/4*KA
  421. 4210  L=0.00508*ZX*12*(LOG(4*12*ZX/DW)-1)
  422. 4220  C=7.36*ZX/(LOG(2*12*ZX/DW)/LOG(10)-0.5)
  423. 4230  KA=(ZX/SQR(L*C)*1000/983.57)
  424. 4240  NEXT W
  425. 4250  Z0=1000*SQR(L/C):XCE=-Z0*TAN(KA*PI/2):CE=-1E+06/2/PI/FMHZ/XCE
  426. 4260  RETURN
  427. 4270  '
  428. 4280  '.....parallel impedance subroutine
  429. 4290  MN=SQR((R7*R8-X7*X8)^2+(X8*R7+X7*R8)^2)
  430. 4300  THETAN=ATN((X7*R8+X8*R7)/(R7*R8-X7*X8))
  431. 4310  IF (R7*R8-X7*X8)<-0.000999999 THEN THETAN=THETAN+PI
  432. 4320  MD=SQR((R7+R8)^2+(X7+X8)^2)
  433. 4330  THETAD=ATN((X7+X8)/(R7+R8))
  434. 4340  THETA=THETAN-THETAD:M=MN/MD
  435. 4350  R9=M*COS(THETA):X9=M*SIN(THETA)
  436. 4360  RETURN
  437. 4370  '
  438. 4380  '.....end
  439. 4390  GOSUB 4430
  440. 4400  GOTO 130   'start
  441. 4410  END
  442. 4420  '
  443. 4430  'HARDCOPY
  444. 4440  GOSUB 4550:LOCATE 25,2:COLOR 14,6
  445. 4450  PRINT " Press 1 to print screen, 2 to print screen & ";
  446. 4460  PRINT "advance paper, or 3 to continue.";:COLOR 7,0
  447. 4470  Z$=INKEY$:IF Z$="3"THEN GOSUB 4550:RETURN
  448. 4480  IF Z$="1"OR Z$="2"THEN GOSUB 4550:GOTO 4500
  449. 4490  GOTO 4470
  450. 4500  FOR QX=1 TO 24:FOR QY=1 TO 80
  451. 4510  LPRINT CHR$(SCREEN(QX,QY));
  452. 4520  NEXT QY:NEXT QX
  453. 4530  IF Z$="2"THEN LPRINT CHR$(12)
  454. 4540  GOTO 4440
  455. 4550  LOCATE 25,1:PRINT STRING$(80,32);:RETURN
  456.